library(tidyverse)
library(sf)
library(leaflet)
## Warning: package 'leaflet' was built under R version 4.0.3
library(htmltools)
## Warning: package 'htmltools' was built under R version 4.0.5
library(htmlwidgets)
library(raster)
library(gstat)
## Warning: package 'gstat' was built under R version 4.0.3
library(spatial)
library(dplyr)
library(jsonlite)
## Warning: package 'jsonlite' was built under R version 4.0.5
library(ggplot2)
library(hrbrthemes)
## Warning: package 'hrbrthemes' was built under R version 4.0.5
library(ggthemes)
options(scipen=999)
##Loading Boundaries
State_Boundaries <- st_read("cb_2018_us_state_500k.kml", quiet = TRUE)
State_Boundaries_GJSON <- st_read("US_State_Boundaries.json", quiet = TRUE)
State_Boundaries_Zip <- st_read("C:/Users/sageg/Desktop/newRrepo/USDA_WLFW_NOBO/StateBoundariesZIP", quiet = T)
##Loading GCT Data
GCT_data <- read_csv("NOBODATA_ForLeaflet_Final.csv")
##Loading RDS
bird_data <- readRDS("NOBO_route_level_trends.rds")
###Adjusting NOBO Trend Datasets
separated_coord_bird_data <- bird_data %>%
mutate(lat= unlist(map(bird_data$geometry,1)),
long= unlist(map(bird_data$geometry,2)))
###Adding Geometries to GCT Data
##Joining Boundaries and GCT Data
GCT_and_Geographies <- st_as_sf(left_join(GCT_data, State_Boundaries_Zip, by =c("STATE" ="NAME")))
## Creating Label
GCT_and_Geographies$label <-
paste("<b>", "<big>", GCT_and_Geographies$STATE,"</b>", "</big>",
"<br>", GCT_and_Geographies$REGION,
"<br>", "Values shown below reflect the combined value of", "<br>", "Original Program Goals and Framework Goals",
"<br>",
"<br>",
"<b> Top 3 Core Conservation Practices: </b>",
"<br>",
GCT_and_Geographies$TOP3CORE,
"<br>",
"<br>",
"<b>", "Financial Assistance:", "</b>","<br>",
"$ ", prettyNum(GCT_and_Geographies$REQFIN_FRAME,big.mark=","),
"<br>",
"<br>", "<b>", "Total CP Coverage, Acres:", "</b>",
"<br>", "Core:", prettyNum(GCT_and_Geographies$ACRE_FRAME_CORE, big.mark = ","),
"<br>", "Supplemental:", prettyNum(GCT_and_Geographies$ACRE_FRAME_SUPP, big.mark = ","),
"<br>", "Core and Supp Combined:", prettyNum(GCT_and_Geographies$ACRE_FRAME_CAS, big.mark = ","),
"<br>",
"<br>", "<b>", "Total CP Coverage, Feet:", "</b>",
"<br>", "Core:", prettyNum(GCT_and_Geographies$FT_FRAME_CORE, big.mark = ","),
"<br>", "Supplemental:", prettyNum(GCT_and_Geographies$FT_FRAME_SUPP, big.mark = ","),
"<br>", "Core and Supp Combined:", GCT_and_Geographies$FT_FRAME_CAS,
"<br>",
"<br>", "<b>", "Total CP Coverage, Number of X:", "</b>",
"<br>",GCT_and_Geographies$X_FRAME,
"<br>",
"<br>", "<b>", "Number of Written Plans:", "</b>",
"<br>",GCT_and_Geographies$WRITTEN_FRAME,
"<br>",
"<br>", "<b>", "Number of Applied Plans:", "</b>",
"<br>",GCT_and_Geographies$APPLIED_FRAME) %>%
lapply(htmltools::HTML)
###Setting up Leaflet with only GCT Data
nobo_summary_leaflet <- leaflet(GCT_and_Geographies)%>%
addProviderTiles(providers$Stamen.Watercolor) %>%
addPolygons(highlightOptions = highlightOptions(fillOpacity = 1, fillColor="cornsilk"),
weight=1,
fillColor= "darkseagreen",
color = "black",
fillOpacity=.5,
label = ~GCT_and_Geographies$label,
labelOptions = labelOptions(
style = list("font-weight" = "normal", padding = "3px 8px"),
textsize = "11px",
direction = "auto")) %>%
addControl("Northern Bobwhite, Grasslands, and Savannas National Partnership 2022-2026: Framework Goal Summary", position = "topright")
## Warning: sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
## Need '+proj=longlat +datum=WGS84'
nobo_summary_leaflet
saveWidget(nobo_summary_leaflet, file="NOBO_SUMMARY_LEAFLET.html")
## Setting up LEaflet
nobo_summary_leaflet <- leaflet()%>%
addProviderTiles(providers$Stamen.Watercolor) %>%
addPolygons(data=GCT_and_Geographies,
highlightOptions = highlightOptions(fillOpacity = 1, fillColor="cornsilk"),
weight=1,
fillColor= "darkseagreen",
color = "black",
fillOpacity=.5,
label = ~GCT_and_Geographies$label,
labelOptions = labelOptions(
style = list("font-weight" = "normal", padding = "3px 8px"),
textsize = "11px",
direction = "auto")) %>%
## Add lci-data
addCircles(data=separated_coord_bird_data,
lng = ~long,
lat = ~lat)%>%
addControl("Northern Bobwhite, Grasslands, and Savannas National Partnership 2022-2026: Framework Goal Summary", position = "topright")
## Warning: sf layer has inconsistent datum (+proj=longlat +datum=NAD83 +no_defs).
## Need '+proj=longlat +datum=WGS84'
nobo_summary_leaflet
#saveWidget(nobo_summary_leaflet, file="NOBO_SUMMARY_LEAFLET.html")
## ggplot
ggplot(bird_data, aes(color=lci_trend)) +
geom_sf()
## leaflet
leaflet(separated_coord_bird_data) %>%
addProviderTiles("CartoDB.Positron") %>%
addCircles(lng = ~long,
lat = ~lat,
color ="red")